Broadband angular spectrum differentiation using dielectric metasurfaces

Signal processing is of critical importance for various science and technology fields. Analog optical processing can provide an effective solution to perform large-scale and real-time data processing, superior to its digital counterparts, which have the disadvantages of low operation speed and large energy consumption. As an important branch of modern optics, Fourier optics exhibits great potential for analog optical image processing, for instance for edge detection. While these operations have been commonly explored to manipulate the spatial content of an image, mathematical operations that act directly over the angular spectrum of an image have not been pursued. Here, we demonstrate manipulation of the angular spectrum of an image, and in particular its differentiation, using dielectric metasurfaces operating across the whole visible spectrum. We experimentally show that this technique can be used to enhance desired portions of the angular spectrum of an image. Our approach can be extended to develop more general angular spectrum analog meta-processors, and may open opportunities for optical analog data processing and biological imaging.

Signal processing is of critical importance for various science and technology fields.Analog optical processing can provide an effective solution to perform large-scale and real-time data processing, superior to its digital counterparts, which have the disadvantages of low operation speed and large energy consumption.As an important branch of modern optics, Fourier optics exhibits great potential for analog optical image processing, for instance for edge detection.While these operations have been commonly explored to manipulate the spatial content of an image, mathematical operations that act directly over the angular spectrum of an image have not been pursued.Here, we demonstrate manipulation of the angular spectrum of an image, and in particular its differentiation, using dielectric metasurfaces operating across the whole visible spectrum.We experimentally show that this technique can be used to enhance desired portions of the angular spectrum of an image.Our approach can be extended to develop more general angular spectrum analog meta-processors, and may open opportunities for optical analog data processing and biological imaging.
Signal processing is vital for several science and engineering disciplines.The increasing demand for high-speed and high-efficiency signal processing has inspired scientists to develop faster, integrated and efficient devices that can process signals and images [1][2][3] .Digital signal processing based on integrated electric circuits can perform complex data processing, but it needs an analog-to-digital converter to discretize the input analog signal into a series of bits.After being processed according to the desired operation, the discrete signals are converted back into analog signals through a digital-to-analog converter.This method suffers from low-speed, large memory requirements and high-power consumption, originating from complicated conversion process 4 and hardware bottleneck 5,6 .
Fourier optics provides a powerful approach to optical processing, by filtering the frequency components of the angular spectrum of the input image.In a conventional 4 F system, the image to be processed is decomposed into its spatial Fourier components by a first optical lens.Various components to filter the spatial frequencies, made by bulky diffractive optical elements, including pinholes, opaque dots, masks and phase plate shifters 57 , can then be exploited for angular spectrum filtering of the incoming waves, leading to imaging applications, such as dark-field 58,59 , schlieren 60,61 , and phase-contrast imaging 62 .However, imparting mathematical operations on the angular spectrum has not been explored to the best of our knowledge.In this paper, we demonstrate metasurfaces that perform mathematical operations on the angular spectrum of an image -"angular spectrum analog meta-processors"-that can be used to perform optical analog processing of the Fourier spectrum-such as differentiation, integration, convolution, and other linear operations as light propagates through them.Depending on the optical operations of interest for angular spectrum processing, we retrieve the required amplitude and phase modulation on the incoming light.We demonstrate differentiation of the angular spectrum based on dielectric metasurfaces across the whole visible frequencies, and we apply it to demonstrate novel applications in the angular spectrum domain, i.e., enhancement of desired angular spectrum features.It should be noted that the angular spectrum differentiation could be also realized by using a lens to extract the angular spectrum, followed by a spatial differentiator to complete angular spectrum differentiation.However, this approach would result in a bulky device.In contrast, here we propose a compact approach to realize angular spectrum differentiation by using a single metasurface.The application potential of enhancing desired portions of the angular spectrum of an image has not been previously exploited in the field of optical analog processing.

Design principles and simulation results
The angular spectrum A p (k x , k y ) of an incoming wavefront can be retrieved by performing the Fourier transform of the in-plane electric where i is the imaginary unit, F indicates the Fourier transform, p = x (y) labels the x-(y-) component, and k x , k y are 2D Fourier domain variables along x and y directions, respectively.A p (k x , k y ) and E p (x, y) represent the complex amplitudes of the electric field in the Fourier and spatial domain, respectively.After passing through an optical diffractive element, the angular spectrum is transformed to ĤA p (k x , k y ), where Ĥ is the angular spectrum operation exerted by the optical diffractive element, modulating the amplitude and phase of A p (k x , k y ) of the incoming wavefront.If an arbitrary order differentiation n j is exerted on the Fourier transform A p (k x , k y ), we have ½C j ðÀixÞ m j ðÀiyÞ n j E p ðx, yÞ expðÀik x x À ik y yÞdxdy ð2Þ where j indicates the differentiation type, and l denotes the maximum number of differentiation type (l is a positive integer, and j ∈ [1, l ]).m j and n j are natural numbers labeling the order of the partial derivative with respect to k x and k y , respectively, and C j is a complex coefficient (associated with the j-th differentiation type), affecting the efficiency and phase delay for the output angular spectrum.Such differential operation in the Fourier domain can be realized by a real-space transfer function t(x, y) = , which is related to Ĥ via Eq.( 2).In other words, we assume that the input image, described by the spatial domain electric field E p (x, y), is filtered by a spatially varying mask with transmission amplitude and phase.We first discuss the general properties of the required mask, and then we show how this real-space transfer function can be achieved with a metasurface.For example, if l = 1, there is only one term in Eq. ( 2), associated with (m j , n j ) = (m 1 , n 1 ) and C 1 .If (m 1 , n 1 ) is chosen as (1, 1), then Ĥ is equal to , associated with t(x, y) = −C 1 xy.The input angular spectrum Ĥ will contain second-order partial derivatives with respect to k x and k y , respectively.If (m 1 , n 1 ) = (1, 0), Ĥ = C 1 ∂ ∂k x with t(x, y) = −iC 1 x, the input angular spectrum will merely consist of the first-order partial derivative with respect to k x .If l = 2, there are two types of differentiation, associated with (m j , n j ) = (m 1 , n 1 ) and (m 2 , n 2 ), and y), a sum over the two firstorder partial derivatives with respect to k x and k y , will be applied on the input angular spectrum.The values of l, C j and (m j , n j ) used for the three types of differentiation can be found in Supplementary Table S1 of Supplementary Note 1.The real-space transfer function t(x, y) provided by the optical diffractive elements should be proportional to xy, x and x + y, respectively, to realize these three types of differentiation.It should be emphasized that, in principle, we can also realize linear combinations of higher-order partial derivatives by flexibly choosing l, C j , and (m j , n j ).For example, we can achieve In the following, we demonstrate metasurfaces tailored to realize these operations.Figure 1a-c schematically show the analog meta-processors that perform the three types of differentiation in the angular spectrum domain.We denote the two opposing edges of the silicon nanopillar cross-section as W u , W v (W u ≥ W v ), namely the long and short sides along u and v directions, respectively.The orientation angle θ represents the angle between the x axis of laboratory coordinate and u axis (Fig. 1d).The local transmission coefficient for ypolarized transmitted light under x-polarized incidence is 63 (see Supplementary Note 2 for the derivation) where t u and t v are the complex transmission coefficients along u and v directions, respectively.Thus, to implement any of the operations discussed above, σ(x, y) needs to match the corresponding real-space transfer function t(x, y).The cross sections of the basic building blocks of each nanopillar can be tailored to control t u and t v in Eq. ( 3), regardless of the operating wavelength.In general, t u and t v are wavelength-dependent, due to the material dispersion and resonant effects within the metasurface unit cells.In general, t u and t v could be chosen to be also position-dependent, in order to increase the flexibility of σ(x, y).Here, to keep the design simple, we have chosen to employ silicon nanopillars with a fixed cross section, so that t u and t v  are wavelength-dependent only.Remarkably, we show that modulating only θ is sufficient to achieve σ(x, y) = t(x, y) for all angular spectrum differential operations considered.Since the required t(x, y) is nonzero at almost every spatial position, the cross section of the silicon nanopillar is chosen to be rectangular, i.e., W u > W v .For a metadifferentiator composed of a square metasurface with a side length a (x ∈ [−a/2, a/2], y ∈ [−a/2, a/2]), we show how to define θ to achieve the aforementioned three types of differentiation.For the first type for the purpose of maximizing the efficiency of the output angular spectrum, corresponding to θ = ±arcsin(4xy/a 2 ) / 2. Using a similar method as for the first meta-differentiator, we can also retrieve θ = ±arcsin(2x / a) / 2 and θ = ±arcsin[(x + y) / a] / 2 for the other two types of meta-differentiators, associated with t(x, y) = −iC 1 x and t(x, y) = −iC 1 (x + y), respectively.The detailed derivation process can be found in Supplementary Note 2. It should be noted that, although the three types of differentiator can work over a wide spectral band, the amplitude and phase shift of the transmissive angular spectrum are wavelength-dependent, as C 1 is dependent on the light wavelength.It should also be noted that the actual meta-differentiator has a limited size, and thus it can only process a certain range of angular spectrum.
Considering an input beam with a Gaussian profile, the waist radius w allowed has an upper limit, which depends on the size of the metadifferentiator.The detailed analysis regarding the image resolution in the Fourier domain can be found in Supplementary Note 3.
The amplitude of the transmissive angular spectrum is proportional to the amplitude difference between t u and t v .Here, the cross-sectional geometry of the silicon nanopillar is optimized to obtain a large amplitude difference between t u and t v at 450, 532 and 685.5 nm.The silicon layer forming the metasurfaces is deposited by magnetron sputtering in the following experiment, and its refractive index and extinction coefficient is measured by an ellipsometer.The optimization of crosssectional geometry can be found in Supplementary Fig. S1 of Supplementary Note 2. The simulated amplitude and phase of σ(x, y) versus θ at four wavelengths are shown in Fig. 1e, f, based on finite difference time domain (FDTD) with the commercial software Lumerical FDTD Solutions.It can be seen that the amplitude of σ(x, y) can be changed from zero to a maximum value as θ is varied from −0.5 π to 0.5 π, and the phase suffers from a π shift as θ changes from a negative value to a positive value.We are able to achieve the aforementioned three types of differentiation by arranging θ so as to fulfill the required distributions of σ(x, y) in a square metasurface.By setting θ equal to arcsin(4xy / a 2 ) / 2, arcsin(2x / a) / 2, and arcsin[(x + y) / a] / 2, respectively, the metasurface distributions and the real parts of σ(x, y)/|σ(x, y)| max required for the three meta-differentiators are shown in Fig. 1g-i, and Fig. 1j-l, respectively.
To verify the angular spectrum differentiation, we use an xpolarized Gaussian beam as the input source to illuminate the metadifferentiators (Fig. 2a), where the associated angular spectrum is shown in Fig. 2b.The corresponding angular spectrum of the filtered beam, given by j ĤA x j 2 =j ĤA x j 2 max , is theoretically calculated for three types of angular spectrum differentiation (Fig. 2c).Three metadifferentiators with 16.8 μm × 16.8 μm in size (e.g., 60 × 60 nanopillars) are used to validate the differentiation functionality.We can first obtain the y-component field intensity distributions after the Gaussian beam passes through the meta-differentiators by FDTD simulations, and then extract their corresponding angular spectrum distributions jA y j 2 =jA y j 2 max by performing the Fourier transform (Fig. 2d).It should be emphasized that the spot waist radius (w = 3.5 μm) is lower than the upper limit of 6 μm.The simulated results are perfectly consistent with the theoretical results, which suggests the image resolution defined in Supplementary Note 3 is rather conservative.Ĥ = C 1 ∂ 2 ∂k x ∂k y results in zero intensity along k y and k x axes, and the maximum intensity is located on the two lines k x + k y and k x − k y (Fig. 2c, d).
results in zero intensity along k y axis and the maximum intensity is located on k x axis only.bandwidth limitations for practically realizing higher order metadifferentiators.The field intensity offered by the meta-differentiators, especially for high order differentiation, is not always consistent with that predicted by theory, e. g., the output field is not symmetric in the long wavelength range, and the field intensity by the meta-differentiators is not so close to zero in the short wavelength range.The detailed study on bandwidth limitation for higher order angular spectrum meta-differentiators can be found in Supplementary Note 5.

Experimental demonstration
We have fabricated three meta-differentiators of 300 μm × 300 μm in size (i.e., 1072 × 1072 nanopillars) on a 500 μm-thick silica substrate to verify the three types of angular spectrum differentiation.Figure 3a-c show scanning electron microscope (SEM) images taken at the lower right corner of each meta-differentiator.The experimental setup (Fig. 3d) is used to check if t(x, y) exerted by the meta-differentiators is consistent with the theoretical values predicted by Eq. ( 2).The polarization directions of two polarizers are mutually orthogonal to extract the cross-polarized transmission field.The distance between objective 1 and the meta-differentiator is adjusted to ensure that the light beam is slightly smaller than the meta-differentiator when reaching the left side of it.Under the illumination of red, green, and blue Gaussian light, the cross-polarized transmission field is recorded by the CCD camera (Fig. 3e-g).The output field intensity distributions are consistent with the theoretical results, as shown in Supplementary Note 6.The field intensity along lines x or y, y and x − y, is zero, and is maximum along lines x + y or x − y, x and x + y, corresponding to the three metadifferentiators with ∂  The meta-differentiator can implement differentiation operation on the angular spectrum of nontrivial objects.The test object is placed between the light source and the first polarizer (Fig. 4a).The distance between the object and objective 1 is first adjusted in order to image the object onto the meta-differentiator.The imaging area is slightly smaller than the meta-differentiator.Then objective 2 is moved to ensure the back focal plane of objective 1 is imaged onto the CCD camera via objective 2. In this scenario, the image recorded by the CCD camera corresponds to the differentiation of the object angular spectrum.The detailed explanation on how to use the experimental setup to retrieve the differentiation of the object angular spectrum can be found in Supplementary Note 7.
In our experiment, the object consists of a stainless steel plate drilled by three parallel rectangular holes and a single circular hole, which has a high transmission contrast between opaque and transparent regions.For the three parallel rectangular holes, the transmitted electric field can be written as The associated angular spectrum and differentiation with the three operations , a single row containing three main spots (Fig. 4c) is transformed into two rows, in which each row contains six main spots (Fig. 4e, f).Instead, when the operation ∂ ∂k x is considered, the three main spots in one row are only transformed into one row of six main spots (Fig. 4g, h).The angular spectrum experiences a first-order differentiation along k x + k y direction with . As a result, under this operation each spot in one row is transformed into two spots along k x + k y direction (Fig. 4i, j).It is worth discussing the antisymmetric profile for the angular spectrum intensity distributions at 450 nm.The polarizer has low extinction ratio of optical powers of perpendicular polarizations at shorter wavelengths, as opposed to longer wavelengths.The y-polarized light beam after polarizer 1 remains and reaches the CCD, which has an intensity comparable to the light beam contributing to angular spectrum differentiation.As the object is a single circular hole, the angular spectrum differentiation operation results are presented in Supplementary Note 9.
Derivative spectroscopy is a technique used to isolate weak spectral features from other unwanted strong features, and it can also be used to strengthen the contrast of a weak band to a strong overlapped signal in the spectral domain.It is widely used for trace analysis, purity test, and stability test 64,65 .While existing techniques operate in real space, our meta-differentiators can be used for isolating features in the angular spectrum domain of an object.If the object angular spectrum, A t (k x , k y ), is much weaker than the additional background angular spectrum, A b (k x , k y ), it is not possible to directly distinguish A t (k x , k y ) and A b (k x , k y ), i.e., one cannot detect the object angular spectrum for trace analysis.However, the object angular spectrum can be detected if the differentiation operation is applied on the object and background angular spectra, and the former differentiation value is much larger than the latter one.The corresponding experimental setup is established as shown in Fig. 5a, where one beam of light illuminates the object (path 1), and then is mixed with another beam of light (path 2), which mimics a strong background.In this case, A t (k x , k y ) is associated with the angular spectrum of the three parallel rectangular holes, and A b (k x , k y ) corresponds to the wide angular spectrum generated by a Gaussian-like source.A t (k x , k y ) can be made much smaller than A b (k x , k y ) by tuning a neutral density filter.In this situation, the mixed angular spectrum without the meta-differentiator around k x = 0 is shown in Fig. 5b, d.It can be seen that the three spots are almost completely invisible (Fig. 5b, d).On the contrary, the six spots that define the angular spectrum decomposition of the target can be instead clearly observed around k x = 0 after the mixed angular spectrum is processed by the meta-differentiator with ∂ ∂k x (Fig. 5c, e).This is because A t (k x , k y ) features a maximum around k x = 0, and the associated differentiation with respect to k x is close to zero, smaller than the differentiation of A b (k x , k y ).
It should be emphasized that our angular spectrum analog metaprocessor is also applicable to higher order differentiation, as has been shown in Supplementary Note 5.Moreover, while here we have focused on differential operations, the meta-processor may be extended to develop many other optical analog processing functionalities in the angular spectrum domain, such as differentiation, integration, convolution and equation solving.Spatial light modulators are also capable of providing spatially varying amplitude and phase distributions, and they are promising for achieving analogous angular spectrum operation functionalities.However, a pixel pitch of spatial light modulators is at least several wavelengths in size, and the resultant angular spectrum processors based on spatial light modulators may generate additional side lobes arising from high order diffraction.In addition to angular spectrum differentiation operation, the metadifferentiators may be explored to realize spatial differentiation operation, when combined with a 4 F system.Based on this functionality, we present high-contrast edge enhancement of resolution test charts and frog egg cells across the whole visible spectrum in Supplementary Note 10.By using nanopillars of varied cross-sections to provide full 0-2 π phase coverage, it is possible to realize other types of angular spectrum differentiation operation, such as It should be emphasized that, while the mathematical operations demonstrated here operate on the complex amplitude of the angular spectrum of the electric field, the CCD camera used for our experiment can only record the field intensity information.Consequently, the phase information cannot be directly retrieved in our experiment.It is also worth discussing how the fabrication error affects the device performance.The silicon film has a surface roughness of 10 nm, and the fabrication uncertainty on the cross section of the silicon nanopillar can be controlled within 10 nm.Our analysis and simulation indicate that the fabrication error hardly influences the spatial resolution and the intensity of the output image.Finally, when assessing the performance of a passive computational metasurface, it is crucial to quantify the throughput efficiency, that is, how the intensity of the output image compares to the intensity of the input image.Different metrics can be used to quantify the efficiency.The "peak efficiency" η peak = max[I out (x, y)] / max[I in (x, y)], is defined as the ratio of the peak intensities in the output and input images, where I out (x, y) [I in (x, y)] is the output (input) intensity map 45,46 .Another metric can be defined by considering the integrated intensities, i.e. η int = RR + 1 À1 I out x,y ð Þdxdy= RR + 1 À1 I in x,y ð Þdxdy.We calculated both metrics for the measurements in Fig. 4. For both metrics, we obtained values in the order of 10 −5 ~10 −3 , depending on the wavelength and the specific operation considered.The relatively low efficiencies are mainly attributed to the cross-polarized configuration of our setup, which leads to a rejection of a large fraction of the input power.The transmission efficiencies can be increased by using highaspect-ratio metasurfaces made by low-loss materials, such as crystalline silicon and TiO 2 .Detailed discussion regarding the transmission efficiency can be found in Supplementary Note 11.

Discussion
In conclusion, we have demonstrated 'angular spectrum analog metaprocessors' that can perform differentiation of the angular spectrum of an object.The optical transfer function can be flexibly designed to target the desired angular spectrum processing functionalities.We have constructed three types of angular spectrum operations based on silicon metasurfaces, and demonstrated their capability of implementing angular spectrum differentiation on practical objects across the entire visible spectrum.We have also presented the use of a metadifferentiator to analyze target objects in the angular spectrum domain.Our results provide a novel approach to image processing in the angular spectrum domain, in contrast to previous strategies in the spatial domain.Such a design not only has broad application prospect for optical analog data processing and biological imaging, but also potentially enables new emerging optical detection techniques.

Fabrication
The sample fabrication process starts with the deposition of 220 nmthick silicon film on a 500 μm silica substrate by using magnetron sputtering.A photoresist layer was spin-coated onto the silicon film and baked in an oven.Afterward, the desired patterns were defined on the photoresist layer by electron-beam lithography (EBL).After development of the photoresist, the pattern is transferred onto the silicon layer by inductively coupled plasma (ICP) etching.The targets with the three parallel rectangle holes and the single circular hole were drilled in the stainless steel plate by etching.

Measurements
The setup to measure the metasurface transmission functions is shown in Fig. 3d.The light from a laser (red: DL-690-010-RS, green: MSL-FN-532-50 mW, blue: DL-450-010-RS) is collected by objective 1 (Olympus UPlanFL N, 10 × /0.30) and the distance between objective 1 and the meta-differentiators is adjusted to ensure that the light beam is slightly smaller than the meta-differentiators when reaching the left side of them.The light field on the meta-differentiators is photographed by objective 2 (Olympus UPlanFL N, 10 × /0.30) and the CCD camera (LBAS-U350-35C).The polarization directions of two polarizers are mutually orthogonal to extract the crosspolarized transmission field.
The setup to measure the angular spectrum differentiation is shown in Fig. 4a.The incident light passes through the object and the object is imaged on the meta-differentiator by objective 1.The object distance and the image distance are adjusted to ensure the image is smaller than the meta-differentiator.Then, the field on the back focal plane of the objective 1 is recorded by objective 2 and the CCD camera.The polarization directions of two polarizers are mutually orthogonal to extract the cross-polarized transmission field.
The setup to measure the angular spectrum isolation is shown in Fig. 5a.The incident light is split into two beams (paths 1 and 2) by the beam splitter 1.The light along path 1 illuminates the object (i.e., the three parallel rectangle holes), and it is then imaged onto the metadifferentiator by objective 1.At the same time, the size of the spot from path 2 is adjusted by a lens (Thorlabs LB1471-A) and the spot is slightly larger than W x .The intensity of the beam from path 1 is attenuated by the neutral density filter to make the corresponding angular spectrum almost completely covered by the wide angular spectrum from path 2. With the same method for measuring the angular spectrum mentioned in Fig. 4, the mixed angular spectrum without and with the metadifferentiators angular spectrum differentiation can be extracted.

Fig. 1 |
Fig. 1 | Angular spectrum meta-differentiator. a-c Angular spectrum analog meta-processors for optical analog differentiation processing in the angular spectrum domain, as light propagates through them.The corresponding angular spectrum operation, Ĥ, and complex real-space transfer function, t(x, y), provided by the silicon metasurfaces are, respectively, / ∂ 2 ∂k x ∂k y and ∝ xy (a), / ∂ ∂k x and ∝ x (b), and / ∂ ∂k x + ∂ ∂k y and ∝ x + y (c).It is assumed that C 1 is equal to C 2 in (c).d Schematic of the rotated silicon nanopillar on a silica substrate with W u = 200 nm, W v = 100 nm, h = 220 nm, and P = 280 nm.e, f Amplitude (e) and phase (f) of σ versus orientation angle θ.The refractive index and extinction coefficient of silicon were measured by an ellipsometer (see Supplementary Fig. S1 of Supplementary Note 2), and the refractive index of silica is extracted from Ref. 66. g-l Schematics of the top view presenting the distribution of the orientation angles (g-i) and distributions of the real parts of the required σ(x, y)/|σ(x, y)| max of the three metadifferentiators (j-l) with ∂ 2 ∂k x ∂k y (g, j), ∂ ∂k x (h, k) and ∂ ∂k x + ∂ ∂k y (i, l), respectively.

Fig. 2 |
Fig. 2 | Angular spectrum intensity distributions for three types of differentiations.a, b Field intensity E x 2 = E x 2 max (a) and angular spectrum intensity A x 2 = A x 2 max (b) of the input x-polarized Gaussian beam with a waist radius of

Fig. 3 |
Fig. 3 | Measured transmission field for three types of differentiations.a-c SEM images of partial samples for angular spectrum differentiation with ∂ 2 ∂k x ∂k y (a) ∂ ∂k x (b) and ∂ ∂k x + ∂ ∂k y (c).The scale bars are 1 μm.d Experimental setup for measurement.The meta-differentiator the CCD camera are on the object and the image planes of objective 2, respectively.e-g Transmission field intensity distributions recorded by CCD with red light (685.5 nm) (e), green light (532 nm) (f) and blue light (450 nm) (g).The scale bars are 100 μm.

∂ 2 ∂k 8 .
The transmitted electric field after the object is plotted in the first panel of Fig. 4b.The associated normalized angular spectrum | A x (k x , k y )| 2 is plotted in the first panel of Fig. 4c, with the intensity profile being shown in the first panel of Fig. 4d.The theoretical angular spectrum differentiations for ∂ 2 ∂k x ∂k y the first panel of Fig.4e, g, i, respectively, where the associated intensity profiles along lines A-B are plotted in the first panel of Fig.4f, h, j, respectively.The angular spectrum at 685.5, 532 and 450 nm are experimentally extracted and recorded by the CCD camera, as shown in the second to the fourth panels in Fig.4b-j.The experiment is well consistent with the above theoretical results.It can be seen that, under the angular spectrum differentiation of ∂ 2 ∂k x ∂k y

Fig. 4 |
Fig. 4 | Measured angular spectrum differentiations of three parallel rectangular holes.a Experimental setup for performing angular spectrum operation of an object.The object and the meta-differentiator are on the object and the image planes of objective 1, respectively.The back focal plane of objective 1 is also the object plane of objective 2 and the CCD camera is on the image plane of the objective 2. b Theoretical (the first panel) and experimental (the second to fourth panels) intensity distributions of three parallel rectangular holes drilled in a stainless steel plate (1 mm thick).The geometrical parameters of the holes are set at W x = 200 μm and W y = 1 mm.c, d Angular spectrum intensity distributions (c) and normalized intensity profiles along lines A-B in (c) (d).e-j Angular spectrum intensity distributions with ∂ 2 ∂k x ∂k y (e), ∂ ∂k x (g), and ∂ ∂k x + ∂ ∂k y (i), and their normalized intensity profiles along lines A-B (f, h, j).

Fig. 5 |
Fig. 5 | Angular spectrum isolation experiment.a Experimental setup for performing angular spectrum isolation at 685.5 nm.b-e Angular spectrum intensity distributions of the mixed light without (b) and with (c) the meta-differentiator with ∂ ∂k x .Their normalized angular spectrum intensity profiles along the lines A-B are shown in (d) and (e), respectively.